function [n, result ] = Euler( funct, x, y, h )
% fungsi untuk mencari solusi PDB dua variabel dengan metode Euler
% fxy = y'
% (x0, y0) --> pasangan nilai awal yang diketahui
% xr --> nilai input untuk tempy yang dicari
% h --> selisih

    x0 = x(1); xr = x(2);
    tempy = y;
    n = (xr-x0) / h;
    
    if (nargin(funct) == 2)
        for r = 1:n
            tempy = tempy + h * funct(x0,tempy);
            x0 = x0 + h;
        end
        result = tempy;
    elseif (nargin(funct) == 3)
        tempy = y(1); tempz = y(2);
        for r = 1:n
            tempz = tempz + h * funct(x0,tempy,tempz);
            tempy = tempy + h * funct(x0,tempy,tempz);
            x0 = x0 + h;
        end
        result = [tempy tempz];
    end
end
